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Dynamics of Vortices in 
Two-Dimensional Magnets 

Franz G. Mertens and A. R. Bishop 



ABSTRACT Theories, simulations and experiments on vortex dynamics 
in quasi- two-dimensional magnetic materials are reviewed. These materials 
can be modelled by the classical two-dimensional anisotropic Heisenberg 
model with XY (easy-plane) symmetry. There are two types of vortices, 
characterized by their polarization (a second topological charge in addi- 
tion to the vorticity): Planar vortices have Newtonian dynamics (even- 
order equations of motion) and exhibit strong discreteness effects, while 
non-planar vortices have non-Newtonian dynamics (odd-order equations of 
motion) and smooth trajectories. These results are obtained by a collective 
variable theory based on a generalized travelling wave ansatz which allows 
a dependence of the vortex shape on velocity, acceleration etc.. An alterna- 
tive approach is also reviewed and compared, namely the coupling of the 
vortex motion to certain quasi-local spinwave modes. 

The influence of thermal fluctuations on single vortices is investigated. Dif- 
ferent types of noise and damping are discussed and implemented into the 
microscopic equations which yields stochastic equations of motion for the 
vortices. The stochastic forces can be explicitly calculated and a vortex 
diffusion constant is defined. The solutions of the stochastic equations are 
compared with Langevin dynamics simulations. Moreover, noise-induced 
transitions between opposite polarizations of a vortex are investigated. 
For temperatures above the Kosterlitz-Thouless vortex-antivortex unbind- 
ing transition, a phenomenological theory, namely the vortex gas approach, 
yields central peaks in the dynamic form factors for the spin correlations. 
Such peaks are observed both in combined Monte Carlo- and Spin Dyna- 
mics-Simulations and in inelastic neutron scattering experiments. However, 
the assumption of ballistic vortex motion appears questionable. 



1 Introduction 



During the past 15 years an increasing interest in two-dimensional (2D) 
magnets has developed. This is a result of (i) the investigation of a wide 
class of well-characterized quasi-2D magnetic materials which allow a de- 
tailed experimental study of their properties (inelastic neutron scatter- 
ing, nuclear magnetic resonance, etc.), and (ii) the availability of high- 
speed computers with the capabilities for simulations on large lattices. 
Examples of these materials are (1) layered magnets like K2CUF4, 
RbaCrCU, (CH3NH3)2CuCl4 and BaM2(X04)2 with M = Co, Ni, . . . and 
X — As, P,...; (2) C0CI2 graphite intercalation compounds Q, (3) mag- 
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netic lipid layers like Mn(Ci8H35 02)2- For the first class of these exam- 
ples the ratio of inter- to intraplane magnetic coupling constants is typically 
10~^ to 10^®. This means that the behavior is nearly two-dimensional as 
concerns the magnetic properties. For the second class, the above ratio of 
the coupling constants can be tuned by choosing the number of intercalated 
graphite layers. For the third class, even monolayers can be produced, using 
the Langmuir-Blodgett method. 

Many of the above materials have an " easy-plane" or XK-symmetry. This 
means that the spins prefer to be oriented in the XY-plane which is defined 
as the plane in which the magnetic ions are situated. The simplest model 
for this symmetry is the 2D classical anisotropic Heisenberg Hamiltonian 
(see Ref. |]) 

H = -J J2 [STS'^ + S;^S^ + il-d)STS^]. (1.1) 

<m,n> 

< m,n > labels nearest neighbors of a 2D lattice; usually a square lattice 
is used. The subscripts x, y and z stand for the components of the classical 
spin vector S; the spin length S is set to unity by the redefinition J — > J/S'^. 
J is the magnetic coupling constant, both ferro- and antiferromagnetic 
materials were investigated. The anisotropy parameter 6 lies in the range 
< (5 < 1; note that (5=1 corresponds to the XY-mode\, not to the so- 
called planar model where the spins arc strictly confined to the XY-plane; 
this confinement is possible only if one is not interested in the dynamics. 



Instead of the exchange anisotropy in the Hamiltonian (1.1), one can also 
use an on-site anisotropy term (S*™)^, which yields similar results. 

There are two kinds of excitations: Spin waves which are solutions of the 
linearized equations of motion, and vortices which are topological collective 
structures. The vortices are responsible for a topological phase transition 
|5| ^ at the Kosterlitz-Thouless transition temperature Tkt- Below Tktj 
vortex-antivortex pairs are thermally excited and destroyed; above Tkt 
these bound pairs dissociate and the density of free vortices increases with 
temperature. It must be noted that an order-disorder phase transition is 
not possible according to the Mermin- Wagner-Theorem 0; the reason is 
that all long-range order is destroyed by the long-wave linear excitations 
in all ID and 2D models with continuous symmetries. 

Both theory and Monte Carlo simulations j|] showed that Tkt is only 
very weakly dependent on the anisotropy 6, except for 6 extremely close to 
0, when Tkt 0. For the materials mentioned above, coupling constants 
were estimated from fits to spin-wave theory, and S values are in the range 
0.01 — 0.6, where Tkt is still close to its value for S = 1. 

However, the out-of-plane structure of the vortices (i. e., the structure 
of the Sz components) depends crucially on S, while the in-plane structure 
{Sx and Sy components) remains the same [|l^. For S > Sc{~ 0.297 for 
a square lattice [Tl}|), static vortices have null Sz components, they are 
termed "in-plane" or planar vortices in the literature. For 5 < Sc there are 
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"out-of-plane" or "non-planar" vortices which exhibit a localized structure 
of the Sz components around the vortex center. This structure falls off 
exponentially with a characteristic length, the vortex core radius [T3| 



in units of the lattice constant a. The core radius increases with diminishing 
5, allowing a continuous crossover to the isotropic Heisenberg model (5 = 0) 
where the topological excitations are merons and instantons fl^ , rather 
than vortices. 

Compared to the vortices in classical fluids Q and superfluids Q, 
there is an important difference: The vortices of the easy-plane Heisenberg 
model carry two topological charges, instead of one: (1) The vorticity q = 
±1, ±2, . . . which is defined in the usual way: The sum of changes of the 
azimuthal angle </> — a,Tcta,n{Sy/ Sx) on an arbitrary closed contour around 
the vortex center yields 27rg; if the center is not inside the contour the 
sum is zero. In the following only the cases q — +1 (vortex) and q ~ — 1 
(antivortex) will be considered, the cases l^l > 1 would be relevant only for 
high temperatures. (2) The polarization or polarity p: For the non-planar 
vortices p = ±1 indicates to which side of the XY-p\a,ne the out-of-plane 
structure points, while p = for the planar vortices (Fig. pi). 
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FIGURE 1. (a) In-plane structure of a static planar vortex (g = -fl). (b) 
Out-of-plane structure of a static non-planar vortex with polarization p = -f 1. 

Both topological charges ^ are crucial for the vortex dynamics, which is 
governed by two different forces: (1) 2D Coulomb-type forces F, which are 



^From the viewpoint of homotopy groups, q and Q — —\qv are tti- and 1x2- 
topological charges, respectively; see Ref. |15|. 
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proportional to the product of the vorticities of two vortices and inversely 
proportional to their distance (assuming that the distance is larger than 
2r^, such that the out-of-plane structures do not overlap), (2) a "gyrocou- 
pling" force |l6[ ^ 

Fg = X X G (1.3) 

where X is the vortex velocity, and X(i:) is the trajectory of the vortex 
center. 



The force (1.3) is formally equivalent to the Magnus force in fluid dy- 
namics and to the Lorentz force on an electric charge. The gyrocoupling 
vector Q 

G = 2nqpe,, (1.4) 

where is the unit vector in z-direction, does not represent an external 
field but is an intrinsic quantity, namely a kind of self-induced magnetic 
field which is produced by the localized S'^-structure and carried along by 
the vortex. For antiferromagnets G vanishes. 

Since the gyrovector ( |l.4| ) contains the product of both topological char- 
ges, the dynamics of the two kinds of vortices is completely different: 

(1) For planar vortices G = 0, therefore they have a Newtonian dynamics 

MX = F, (1.5) 

where the vortex mass M will be defined in section 2.2. However, in the 
simulations the trajectories are not smooth due to strong discreteness ef- 
fects 

(2) Non-planar vortices have smooth trajectories if the diameter 2r„ of 
the out-of-plane structure is considerably larger than the lattice constant 
[TqI ; i. e., if (5 is not close to Sc- For steady state motion, the dynamics is 
described by the Thiele equation [^6[ |l^ 

XxG = F, (1.6) 

which was derived from the Landau-Lifshitz equation (section 2) for the 
spin vector S™(t). This equation is identical to Hamilton's equations with 
Hamiltonian (1.1). 

However, for arbitrary motion the Thiele equation (1.5) is only an ap- 
proximation because a rigid vortex shape was assumed in the derivation. 
If a velocity dependence of the shape is allowed pfl, ElJ , an inertial term 



AfX appears on the 1. h. s. of Eq. (1.6). However, this 2' -order equation 



could not be confirmed by computer simulations ||2^, g4|. The reason will 
be discussed in section 2; in fact, a collective variable theory reveals that 
the dynamics of non-planar vortices can only be described by odd-order 
equations of motion. Therefore the dynamics is non- Newtonian. 

Another interesting topic is the dynamics under the influence of thermal 
fluctuations (section 3). Here white noise and damping are implemented in 
the Landau-Lifshitz equation. The same type of collective variable theory 
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as in section 2 then yields stochastic equations of motion for the vortices. 
The stochastic forces on the vortices can be calculated and the solutions of 
the equations of motion are compared with Langevin dynamics simulations. 
In this way, the diffusive vortex motion can be well understood. 

For somewhat higher temperatures (still below Tkt) the thermal noise 
can induce transitions of non-planar vortices from one polarization to the 
opposite one. Theoretical estimates of the t rans ition rate are compared 
with Langevin dynamics simulations (section 3.3). 

For the temperature range above Tkt, so far only a phenomenological 
theory exists, the vortex-gas approach (section ^). It contains only two 
free parameters (the density of free vortices and their r.m.s. velocity). The 
theory predicts "central peaks" in the dynamic form factors for both in- 
plane and out-of-plane spin correlations. Such peaks are observed both 
in computer simulations and in inelastic neutron scattering experiments, 
and many of the predicted features are confirmed. However, a recent direct 
observation of the vortex motion in Monte Carlo simulations has revealed 
that a basic assumption of the theory, namely a ballistic vortex motion, 
might not be valid. (In section || we confine ourselves to the ferromagnetic 
case, although easy-plane antiferromagnets were also investigated). 

The conclusion in section ^ will discuss the ingredients of a theory which 
can fully explain all observations above Tkt- 



2 Collective Variable Theories at Zero 
Temperature 

2.1 Thiele Equation 

The spin dynamics is given by the Landau-Lifshitz equation, which is the 
classical limit of the Heisenberg equation of motion for the spin operator 

^--S'-x^ (17) 
where H is the Hamiltonian, in our case that of the anisotropic Heisenberg 



model (LI). The meaning of Eq. (L7) is that the spin vector S™ pre- 



cesses in a local magnetic field B, with cartesian components Ba = — 
In spin dynamics simulations (see below), the Landau-Lifshitz equation is 
integrated numerically for a large square lattice, typically with 72 x 72 lat- 
tice points. As initial condition a vortex is placed on the lattice and the 
trajectory X(t) of the vortex center is monitored. 

This trajectory is then compared with theory, i. e. with the solution of an 
equation of motion for the vortex. The standard procedure to obtain this 
equation consists in taking the continuum limit S™(i) — S(r, t) where r is 
a vector in the XY plane, and in developing a collective variable theory. 
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The simplest version makes the travelhng wave ansatz |l^ 

S(r,<) = S(r-X(t)), (1.8) 

where the functions Sa on the r. h. s. describe the vortex shape. (Strictly 
speaking, these functions should bear an index to distinguish them from 
the functions on the 1. h. s.) 

As the equation of motion is expected to contain a force, the following 
operations are performed with Eq. (1.7) which yield force densities 

as 



dXi dt 



= -S 



SH dS 



SS dX, 



dX, 
(1.9) 



with i ~ 1,2 and Hamiltonian density Ti,. According to the ansatz (1. 



dS 

dt 



dXf ' 



(1.10) 



is inserted, with summation over repeated indices. Integration over r then 
yields the equation of motion 



with the force 



GX = F, 



and the gyrocoupling tensor 



G,;, = / d^r S 



9S 9S 



dX^ dXi 



= d 



dX,, ' 



d(f) dip dcj) dip 



dX^ dXj dXj dX^ 



(1.11) 
(1.12) 

(1.13) 



The expression on the right was obtained by introducing the two canonical 
fields 

= arctan(S';^/S':r), "0 = & (1.14) 

instead of the three spin components Sa with the constraint 5' = 1. The 
static vortex structure is ||l^ 



q arctan — 



1-at 



for r' <^ r„ 



ipQ = pa2\ p-Te for r' » r„ 



where the constants oi and 02 can be used for a matching, and 

r' = r - X . 



(1.15) 
(1.16) 
(1.17) 

(1.18) 



1. Dynamics of Vortices in Two-Dimensional Magnets 7 



The integral then yields 

Gy = Ge^J , G = 2TTqp (1.19) 

where is the antisymmetric tensor. Interestingly, only the value p of the 
Sz component at the vortex center enters the final result [p^ , i. e., the 
out-of-plane vortex structure needs not to be explicitly known here. 

As G is antisymmetric, Eq. ( |l.ll[ ) is identical to the Thiele equation ( |l.6| ) 
with gyrovector G = Ge^. Since p = and thus G = for planar vortices, 
the Thiele equation is incomplete in this case. Obviously there must be a 



non- vanishing term on the 1. h. s., this will be an inertial term (section 2.2) 



Moreover, the ansatz (1.8), and thus the Thiele equation, is only valid for 
steady-state motion when the vortex shape is rigid (in the moving frame). 
This includes not only translational motion with constant velocity Vq but 
also rotational motion with constant angular velocity luq. Both types of mo- 
tion can be obtained by considering two non-planar vortices at a distance 
2i?o which drive each other by their Coulomb force F = 2nqiq2/ {2Rq). 
For a certain velocity this force is compensated by the gyrocoupling force 
(^). In fact, since each vortex carries two charges, there are four physi- 
cally different scenarios which represent stationary solutions of the Thiele 
equation. They fall into two classes: If the gyrovectors of the two vortices 
are parallel (i. e., qipi = q2P2), a vortex-vortex or vortex-antivortex pair 
rotates with luq on a circle with radius Rq around each other, where 

2 pi 2 p2 

If the gyrovectors are antiparallel, the pair performs a parallel transla- 
tional motion with velocity Vq and distance 2i?o, where VqRo = q2/{2pi) = 

-gi/(2p2). 

Both kinds of scenarios also appear for vortices in (super)fluids |l^, 
However, in these contexts there are only two physically different situations: 
vortex-vortex rotation and vortex-antivortex translation. 



2.2 Vortex Mass 

The above scenarios for magnetic vortices were not tested by computer 
simulations until 1994, and the results were very surprising ||2^ , |2^ : Using 
a square or circular system with free boundary conditions, two of the four 
scenarios (vortex-vortex rotation and vortex-antivortex translation) were 
confirmed, but not the two other ones: For vortex-antivortex rotation and 
vortex-vortex translation the observed trajectories showed pronounced os- 



cillations around the trajectories predicted by the Thiele equation (1.6) or 



( ll.llD , see Fig. g. 

In fact, such oscillations had already been predicted j2^, |l| by assuming 



that the vortex shape in the travelling wave ansatz (1.5 ) depends explicitely 
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FIGURE 2. Vortex-antivortex rotation on a circular system with radius L = 36 
and free boundary conditions. Only the trajectory of the vortex center is plotted, 
the center of the antivortex is always opposite. The mean trajectory is a circle 
with radius Ro = 15.74. 



on the velocity X. This leads to an additional term -^^^j in Eq. ( 1.10 ) 
and thus to an inertial term in the Thiele equation 

MX + GX = F , 

where M is the mass tensor 

dS /" \ d4> dip d4> dip 1 



(1.21) 



d^.S^x 



dX, dX, dX, dX, 



(1.22) 



This integral can be easily evaluated if the vortex is placed at the center of a 
circular system with radius L and free boundary conditions. The dominant 
contribution stems from the outer region Uc < r < L, with Uc — 3r„, where 
the velocity dependent parts of the vortex structure are |]lO[ |^ 

q 



^1 



;ix2Xi 



~XiX2) 
X2X2) 



(1.23) 
(1.24) 



01 = p{xiXi 

Together with the static parts (pQ from Eq. (1.15) and V'o, which falls off 
exponentially outside the core (pT^), the rest mass M is obtained 

"2 L 



M,,. 



M 



4(5 ttr 



+ Cm. 



(1.25) 
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The constant Cm stems from the inner region < r < Oc, where the 
vortex structure is not well known due to the discreteness, but Cm is not 
important for discussion. There is also a velocity dependent contribution to 
the mass which is negligible because the vortex velocities in the simulations 
are always much smaller than the spin wave velocity, which is the only 
characteristic velocity of the system. 

The above vortex mass is consistent with other results in the literature: 
Generalizing the momentum of solitons in ID magnets ||25|l , the vortex 
momentum 

P = ~ j A^rtjjVcj) (1.26) 

is defined and can be shown to be a generator of translations j2l). Then 
P = MX results, but P is not a canonical momentum because the Poisson 
bracket {Pi,P2} = C does not vanish. For the kinetic energy AfX^/2 is 
obtained [0, therefore this energy and the rest energy E = tt ln{L/ac)+CE 
both show the same logarithmic dependence on the system size L. 

For a continuum model of a 2D ferromagnet with uniaxial symmetry and 
a magnetostatic field, a slightly different vortex momentum was defined 
[ p6| , namely 

PPT= /d^rrxg, (1.27) 



where g = V<^ x Vip represents the gyrovector densit y, cf. Eg. (1.13) which 



i s rela ted to the gyrovector as described below Eq. ( 1.1£ ). The definition 



( 1.27 ) depends on the choice of origin of the system and is not a generator 



of translations. Nevertheless, if the time derivative of PpT is set equal to 



the force F on the vortex, the generalized Thiele equation ( 1.21 ) is again 
obtained Therefore the vortex dynamics is qualitatively the same, as 
will be discussed now. 



Eq. ( 1.21 ) has the same form as that for an electric charge e in a plane 
with a perpendicular magnetic field B and an in-plane electric field E. I. e., 
the vortex motion is completely analogous to the Hall effect: The trajectory 
is a cycloid with frequency 

(1.28) 

cf. the cyclotron frequency eB/{Mc), where c is the speed of light. It is 
possible to transform to a frame where the force F vanishes and the vortex 
rotates (i. e., guiding center coordinates ^). 

The cycloidal trajectories can explain qualitatively the oscillations ob- 
served in the simulations (Fig. ^). Moreover, the fact that oscillations are 
observed only for two of the four scenarios for two vortices (see above) can 



also be explained: A generalization of the ansatz (1.8) to the case of two 
vortices yields two coupled equations of motion with a 2- vortex mass tensor 
[ p8| p3| . In the case of vortex- vortex rotation and vortex-antivortex trans- 
lation there are cancellations in this mass tensor which prevent oscillations, 
in agreement with the simulations. 
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However, a quantitative comparison of Eq. 
reveals two severe discrepancies 



(1.21) and the simulations 



24 



(1) The mass M = G/lOc which is obtained from Eq. ( |1.28 ) by inserting 
the observed frequency, turns out to be much larger than predicted by 
Eq. (1.25). Moreover, the dependence on the system size L is linear, while 



( 

Eq. (1.25) predicts a logarithmic dependen ce. 

(2) Instead of the one frequency Uc of Eq. (|l.28| ), two frequencies ijJi^2 are 
observed in the spectrum of the oscillations (Fig. As wi and ui2 are close 
to each other, a pronounced beat is observed in the trajectories (Fig. ^): 
The cycloidal frequency is (a-'i+a-'2)/2, but the shape of the cycloids changes 
slowly with [uj2 — wi)/2. 
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FIGURE 3. Upper panel: Fourier spectrum of the radial displacements 
r{t) = R{t) — i?o of the vortex in Fig. |2| from the mean trajectory; the inte- 
gration time is 20000 in units (J5)~^. Lower panel: Spectrum of the azimuthal 
displacements <p(t) = 4>(.t) —wot from the mean trajectory, where wo is the angular 
velocity. 

Both discrepancies have nothing to do with 2-vortex effects, because they 
also occur in 1-vortex simulations . Here the vortex is driven by image 
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forces, in analogy to electrostatics. The most convenient geometry is a 
circular system with radius L and free boundaries. In this case there is only 
one image vortex which has opposite vorticity, but the same polarization 
[ p2| [2g| ]. The radial coordinate of the image is L^/i?, where R is the vortex 
coordinate. The vortex trajectory can be fitted very well to a superposition 
of two cycloids. This observation will be explained in the next section. 



2.3 Hierarchy of equations of motion 



Recently a very general collective variable theory was developed for non- 
linear coherent excitations in classical systems with arbitrary Hamiltonians 
[ p^ . In this theory the dynamics of a single excitation is governed by a 
hierarchy of equations of motion for the excitation center X(t). The type of 
the excitation determines on which levels the hierarchy can be truncated 
consistently: "Gyrotropic" excitations are governed by odd-order equations 
and thus do not have Newtonian dynamics, e. g. Galileo's law is not valid. 
"Non-gyrotropic" excitations are so-to-speak the normal case, because they 
are described by even-order equations, i. e. by Newton's equation in the first 
approximation. 

Examples of the latter class are kinks in ID models like the nonlinear 
Klein- Gordo n fa mily and the planar vortices of the 2D anisotropic Hcisen- 
berg model ( IT ) . The non-planar vortices of this model represent the sim- 
plest gyrotropic example. 3D models have not been considered so far. 

The basis fo the above collective variable theory is a ''''generalized trav- 
elling wave ansat^^ for the canonical fields in the Hamiltonian. For a spin 
system, as considered in this review, this ansatz reads 



S(r,t) = S(r-X,X,X,...,X(")). 



(1.29) 



This generalization of the standard travelling wave ansatz (1.8) yields an 
{n + l)'^-order equation of motion for X(t), because Eq. (1.10) is replaced 

by 



dS 

dt 



OS 

dX, 



ds_ 

dX, 



-X, 



dX 



'9S (n+l) 
(n) 



(1.30) 



3 



The same procedure as described below Eq. ( 1.10| ) then yields the {n+l)*^- 
order equation. 

The case n — 1 corresponds to the 2"'^-order equation (1.21). In one 



dimension the antisymmetric tensor G from Eq. (1.13) vanishes and (1.21) 
reduces to a Newtonian equation. The same is true for the 2D planar vor- 
tices where dj = because p = in Eq. (1.19). 
The case n — 2 yields the 3'''*-order equation 



AX + MX + GX = F(X) 



(1.31) 
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with the 3"^ -order gyrotensor 



A, . / d^.sH X ^ . U^r - ) . (1.32) 



For the above ID models and the 2D planar vortices nothing changes be- 
cause Aij turns out to be zero (below). But the 2D non- planar vortices are 
the first gyrotropic example. 



For the calculation of the integral (1.32) the acceleration dependence of 



the outer region of the vortex is needed |24 



V'2 = ^ (xiXi + a;2X2) (1.33) 

02 = |j(x2li-xiX2) In^- (1-34) 
ad eL 



Together with the static parts (po and ipo this yields 

_G 
16"^ 



Ay = Aey , —{L^ -al) +Ca, (1.35) 



where the constant Ca stems from the inner region. 

In the simulations the dynamic in-plane structures 0i and 02 cannot be 
clearly observed because it is difficult to subtract the static structure (pQ 
which varies drastically with the vortex position. However, the dynamic 
out-of-plane structures ■01 and ^02 can be observed and distinguished by 
looking at specific points of the trajectory: E. g., at the turning points in 
Fig. ^ the acceleration is maximal while the velocity is small. Figs. |5[ |6| 
confirm the structure of 0^2 in the outer region, they also show that this 
structure oscillates as a whole, with frequency (lui + ix'2)/2. An alternative 
interpretation of this oscillating structure in terms of certain "quasilocal" 
spinwaves will be given in section 2.4. 

As 02 and 02 depend linearly on Xj and as Aij contains derivatives with 
respect to Xj , the leading contribution to Aij is independent of velocity and 
acceleration. There are higher-order terms which are negligible, though; cf. 
the discussion of the mass. For these reasons the 1. h. s. of the equation of 



motion ( 1.31 ) is linear. 

The radial Coulomb force F on the r. h. s. can be linerarized by expanding 
in the small radial displacement r(t) — R(t) — Rq from the mean trajectory 
which is a circle of radius Rq 

F{R)^F{Ra) + F'iRo)r. (1.36) 



The 3 -order equation of motion (1.31) then has the following solutions 



24 : A stationary solution and a superposition of two cycloids 

r(t) = aiCOSLUit + a2COSLU2t (1-37) 
RaLp(t) = 6i sinwii 4- 62 sina;2t 
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5 10 15 20 25 

FIGURE 4. First part of the trajectory of a vortex with g = p = 1 on a circular 
system with radius L = 36 and free boundary conditions. The small diamonds 
(o) mark the position of the vortex in time intervals of 10(JS)~^. 



where if = (j) — ujot is the azimuthal displacement. The results for wi,2 can 
be considerably simplified for i?o ^ which is the case in the simulations. 
The frequencies loi,2 form a weakly-split doublet. The mean frequency de- 
pends on A, but not on M: 



Lo — ^ui\tjj2 = \/ G/A ^ 1/L 
for large L. Contrary to this, the splitting of the doublet 



Aw ~ L02 — UJl = M / A 



~1X 



(1.38) 



(1.39) 



is proportional to the mass. 

Ks G — 2TTqp is known from Eq. (1.1£), the last two equations are suf- 
ficient to determine M and A by using the frequencies ^ observed in 
the simulations. In Ref. the data for lattice sizes L = 24 ... 72 were 
extrapolated for _Ro ~^ by using several runs for small Rq/L, because the 
theoretical predictions (1.25) and (1.35) were made for Rq = 0. For the 
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FIGURE 5. Out-of-plane structure of the vortex at the turning point of the 
trajectory in Fig. ^ Here the acceleration has a maximum and points in the radial 
direction, while the velocity is small and points in the azimuthal direction. 

anisotropy (5 = 0.1 the data for A are well represented hy A = CL" + Aq 
with a = 2.002, C = 4.67 and = 40. This agrees rather well with the 



i^-term in (1.35), the constant Ca cannot be tested. However, M « 15 
is nearly independent of L, in contrast to the logarithmic dependence in 
Eq. ( 1.25| ). This can be connected to the observation that the velocity 



dependent part V'l seems to approach an L-independent constant at the 



boundary, in contrast to the predicted l/r-decay in Eq. (1.23) 



A very interesting point is the discussion of the size dependence of the 



different terms in the equation of motion (1.31). As every time derivative 
of Xi contributes a factor of 0Ji^2 — 0{l/L), the orders of the terms are 

AX^ = 0{1/L),MX^ = 0{\nL/L^) and GX, = 0{l/L). Therefore the 
strong 3'''^-order term cannot be neglected when the weak 2"'^-order term 
is retained. This is the reason for the two severe discrepancies between 



the simulations and the predictions of the 2"'^-order equation ( 1.21 ), which 
were discussed in section 2.2. However, the neglection of both the 2"^^- and 
3'^'^-order terms represents a consistent first approximation, namely, the 



Thiele equation (1.11) 



The next consistent approximation in the hierarchy is the 3"^-order equa- 
tion ( 1.3l|), as discussed above. An investigation of even higher-order terms 



in Ref. ||2J| confirms the conjecture that only the odd-order equations of 
the hierarchy represent valid approximations for gyrotropic excitations (For 
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FIGURE 6. Same as in Fig. |5|, but at the 8* turning point, where the accelaration 
points in the negative radial direction. 



non-gyrotropic ones all members of the hierarchy are even). 

The odd higher-order equations predict additional frequency doublets 
1^3,41 <^5.6 etc. These frequencies normally vanish in the background of the 
spectrum. However, they become visible if the simulations are specially 
designed: For a vortex in the center of a quadratic system with antiperiodic 
boundary conditions, all image forces cancel exactly and only the small 
pinning forces remain. For this configuration two additional doublets with 
strongly decaying amplitudes were indeed observed [p^ . 

Finally we mention that cycloidal vortex ti^aiectories have been found 
not only in field-theoretic models for magnets but also for (1) neu- 

tral and (2) charged superfiuids: 

(1) In the Ginzburg-Landau theory, which describes vortex motion in thin, 
neutral, superfiuid films, the usual assumption of incompressibility was 



dropped |31 . A moving vortex then exhibits a density profile in the region 
outside the vortex core. This profile is similar to the velocity dependent out- 
of-plane structure (1.23) of the magnetic vortices, and the consequences are 
also similar: There are small-amplitude cycloids |32| superimposed on the 
trajectories in both 2-vortex scenarios (vortex-vortex rotation and vortex- 
antivortex translation, see end of section pTl| ). 

(2) For the dynamics of vortices in a charged superfiuid the same kind of 
superimposed cycloids was found, again for both 2-vortex scenarios | |3^ . 
Here a field-theoretic model was used, where a charged scalar field is min- 
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imally coupled to an electromagnetic field and a 0*-potential is included; 
this is proposed as a phenomenological model for a superconductor. 

However, compared to the dynamics of magnetic vortices as reviewed in 
this article, there is an important difference: For all the models discussed 
above, there is only one cycloidal frequency; correspondingly the vortex 
dynamics is governed by 2"'^-order equations of motion. 

It will be interesting to study many other physical contexts in which 
topological vortex structures appear (e. g., dislocations in solids and flux 
lattices, vortex filaments in compressable fluids and complex Ginzburg- 
Landau models). We expect that the details of the vortex dynamics differ, 
depending on the order of the equations of motion, for instance (see begin- 
ning of this section). 



2.4 An alternative approach: Coupling to magnons 



The spirit of the generalized travelling wave ansatz (1.29) differs consider 



ably from the well-known standard ansatz, which reads for a spin system 
S(r,t) = So(r-X(t))+x(r,t). (1.40) 

Here So represents the static structure of a single nonlinear coherent exci- 
tation (a vortex in our case) and x represents a magnon field (or meson 
field for other kinds of systems). Since So is static, but the shape of the 
excitation usually depends on the velocity, a part of the dynamics obviously 
must be taken over by the spin waves. 

Based on earlier work, this concept has very recently been carried through 



40 . First, the magnon modes in the presence of a single static vortex were 



obtained by a numerical diagonalization for relatively small, discrete sys- 



tems with fixed (Dirichlet) boundary conditions (BC) 35, S^]. Analyti- 
cal investigations were done for planar vortices in antiferromagnets ^ 
and ferromagnets ||3^, These articles demonstrated nontrivial proper- 
ties of the eigenmodcs, c. g., the presence of quasi- local (resonance type) 



34 , |35| or truly local |36 modes. Moreover, the relevance of these modes 
for the vortex dynamics was shown, in particular the transition between 
planar and non-planar vortices was investigated ||3j, |3^ . 

However, for non-planar vortices all this was based on numerical diago- 
nalization. But very recently, a general theory of vortex-magnon coupling 
was developed for arbitrarily large systems with circular symmetry and 
general BC The S-matrix for vortex-magnon scattering was calculated 
and expressed by Bessel and Neumann functions. Using the S-matrix, gen- 
eral formulas for the eigenfrequencies were obtained, as a function of the 
parameters and size of the system, and for different BC, namely for Dirich- 
let (fixed) and von Neumann (free) BC. 

There is a very good agreement between the frequencies of the two lowest 
quasi-local modes and the frequencies uji,2^ which were observed for the 
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vortex trajectories in spin dynamics simulations for the discrete system 
(section (2.2). The error is only 0.8 % for the mean frequency Co = y^oJjll^ 



and 4 % for the splitting Acj = L02 — ^i- This demonstrates very clearly 
that the cycloidal vortex motion is accompanied by certain magnon modes, 
namely by quasi-local modes. These modes are both extended (like the other 
magnon modes) and local, because they exhibit a localized structure around 
the vortex center (Fig. 0). The structure of these modes is qualitatively 
similar to the dynamic vortex structure, which oscillates with the rythm 
of the cycloidal motion (Figs. ||, ^). However, only the out-of-plane part tp 
of this stru cture can be clearly observed, and only in the outer region (see 
section 2.3). 







FIGURE 7. One of the two low-lying quasi-local spin wave modes in a complex 
representation, obtained by a numerical diagonalization for a circular system with 
a vortex in its center, from (In this reference m corresponds to our tp). 

As the frequencies lui^2 in the cycloidal vortex motion are the same as 
the frequencies of the lowest quasi-local modes, and as the latter ones were 
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calculated analytically, the parameters A and M in the 3'''^-order equation 
of motion ( I.31| ) can also be calculated analytically |^ . The results can be 
expressed by the first root xi of the equation aJi{x) + bxJ[{x)r,j/L = 0, 
where Ji is a Bessel function. The values for a and b define the boundary 
conditions: a = 0, 6 = 1 for free BC and vice versa for fixed EC. 
The general result 

in the limit of large L specializes to A ~ 4.634i^ for 5 = OA and free BC, 
which agrees perfectly with the result 4.67L^ '^"^ from the simulations (sec- 
tion 2.3). For fixed BC A — 1.07i^ is about four times smaller. T he ge neral 
result for the vortex mass M is much more complicated than Eq. ( 1.41 ) and 
can be found in Ref. The numerical value M = 14.74 for S = 0.1 and 
free BC agrees well with the value 15 from the simulations (section 2.3). 
However, both results are independent of the system size L, in contrast to 
the logarithmic L-dependence ( 1.25 ) predicted by the integral ( 1.22 ) in the 
collective variable theory, cf. the discussion in section 2.3. For fixed BC, 
AI — 7.661 is obtained, about half of the above value for free BC. 

The fact that both the mass M and the factor A in the 3"^-order gyro- 
coupling term depend strongly on the boundary conditions appears very 
natural because the vortices are not localized excitations like solitons in 
ID; the in-plane vortex structure falls off with l/r. 

It is important to note that the higher-lying quasi-local modes could be 
calculated by the same methods, and are expected to agree with the higher- 
order doublets uj3^4,uj5fi, ■ ■ ■ appearing in the spectrum of the trajectories. 
Each additional doublet corresponds to taking into account two additional 
orders in the generalized travelling wave ansatz (1.29), as discussed in sec- 



tion 2.3 



Finally, we mention the recent development of a very general collec- 
tive variable theory for an arbitrary Hamiltonian H[(j3,^] supporting non- 
linear coherent excitations, without making any approximations ||2^, 
This theory starts with the standard ansatz for 4>{r,t) and 'i/;(r,<), cf. 
Eq. ( 1.40| ), but imposes constraints between the meson fields and the func- 
tions (j)o{r-,Xi, . . . , Xm) and ipo{r; Xi, . . . , Xm), which describe the shape 
of the excitation. Here the Xi{t) are m collective variables for the position 
and the internal modes of the excitation (if there are any). There are then 
2m constraints in order to preserve the correct number of degrees of free- 
dom. Consequently, the mathematical formalism is based on the classical 
limit of Dirac's quantum mechanics for constraint systems. The equations 
of motion for the collective variables are a generalization of Thiele's equa- 
tion (1.11). The rank of the gyrocoupling matrix G is used to classify the 
excitations: In the case of vanishing G the excitations have an effective 
mass and exhibit Newtonian dynamics; in the case of regular G, the ex- 
citations behave like charged, massless particles in an external magnetic 
field, similar to the gyrotropic excitations defined in section 2.3. The above 
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theory is a generalization of earlier work of Toniboulis [Q, Boesch et al. 

and Willis et al. |Q, which applies only to "standard" Hamiltonians 
(i. e., consisting of the sum of kinetic and potential energy terms). This 
generalization is necessary, e. g., in the case of magnetic systems which 
cannot be modelled by standard Hamiltonians. 



3 Effects of thermal noise on vortex dynamics 

3.1 Equilibrium and non- equilibrium situations 

As already mentioned in the introduction, below Tkt vortex-antivortex 
pairs appear and vanish spontaneously due to thermal fluctuations. But 
these pairs do not move and therefore give no contribution to the dynamic 
correlation functions. Above Tkt, some of the pairs unbind and the free 
vortices can move. In section 4, this situation will be investigated by a 
phenomenological theory, namely the vortex gas approach. 

On the other hand, the effect of thermal fluctuations on single vor- 
tices can be studied by putting a vortex in a thermal bath. This is a 
non-equilibrium situation, in fact it takes a long time until equilibrium 
is reached: The vortex very slowly approaches the boundary where it anni- 
hilates together with an image antivortex; during this process spin waves 
are radiated which arc eventually thermalized. 

The vortex motion with thermal noise is a random walk process, where 
a vortex diffusion constant can be defined. This offers the possibility to 
develop an ab-initio theory for the dynamic form factors. This can be 
compared with the phenomenological vortex gas approach, which assumes, 
however, a ballistic motion (section^). 

3.2 Collective variable theory and Langevin dynamics 
simulations 

In principle, the generalization of a collective variable theory to finite tem- 
peratures is a straightforward procedure consisting of 4 steps; 

(1) Thermal noise and damping (because of the fluctuation-dissipation the- 
orem) are introduced into the microscopic equations. 

(2) A travelling wave ansatz for a nonlinear coherent excitation is made 
which yields equations of motion with stochastic forces acting on the exci- 
tation. 

(3) The solutions of these stochastic o. d. e.'s are used to calculate the 
variances of the trajectory, which contain as a factor an effective diffusion 
constant for this excitation. The dependences of this constant on the tem- 
perature and other parameters can be discussed. 

(4) The predicted variances and the effective diffusion constant are com- 
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pared with the same quantities as observed in Langevin dynamics simula- 
tions. 

Step 1: Introduction of thermal noise 

This step is very problematical, although there are many papers in which 
the problems are either not appreciated or hidden. There are basically two 
major problems: 

(a) The microscopic equations can often be written in different ways which 
are equivalent. In the case of spin systems, the Landau-Lifshitz equation for 
the spin S is equivalent to the Hamilton equations for (p and ^p. However, 
if noise is implemented, e. g., by an additive term, the results may be 
qualitatively different (see below). 

(b) Typically, either additive or multiplicative noise can be used, but many 
papers do not give a reason why one of the two types was chosen. However, 
the results are usually qualitatively different for the two types of noise. 

In the case of the vortex dynamics, both major problems have been 
investigated in a preprint p5[| : Additive noise in the Hamilton equations 
yields a vortex diffusion constant Dy which depends logarithmically on 
the system size, while in the multiplicative case Dy is independent of the 
system size. Additive noise in the Landau-Lifshitz equation ^ yields a 
diffusion constant with the same logarithmic term as in the case of additive 
noise in the Hamilton equations, but the small constants that appear in 
addition to the logarithmic term differ for the two cases because the vortex 
core gives different contributions. 

Unfortunately, the additive noise in the Landau-Lifshitz equation has 
an unphysical feature, namely the spin length is not conserved. In the 
collective variable theory this problem is overcome by a constraint, while 
in the simulations a renormalization of the spin length is necessary after 
every time step|^. 

However, taking additive noise and using these little tricks is not really 
necessary, because there is a better noise term which eventually leads to the 
same results, but which is well motivated on physical grounds and which 



conserves the spin length In the Landau-Lifshitz equation (1.7) 

each spin S™ precesses in a local magnetic field B which is the gradient 
of the energy with respect to the spin compononts. This local field is the 
only way through which the spin S™ can feel any changes in its environ- 
ment. Adding a thermal noise term to the local field thus accounts for the 
interaction of the spin with magnons, phonons and any other thermally 
generated excitations. The stochastic Landau-Lifshitz equation then reads 



^Technically, this is achieved by adding a Lagrange parameter multiplying the 
constraint |^6|. This means that a multiplicative noise term appears besides the 
additive one. 
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dS^ 



' dH 



h™(t) 



dS'" 



(1.42) 



Because of the cross product with S™ this noise is multiplicative. In Ref. 
[|5| Gaussian white noise is used: 



< h^m;,{t') > 



> = 



D5"''^6a,3S{t - 1') 



(1.43) 
(1.44) 



where D = 2ekBT is the diffusion constant and a, (3 denote cartesian com- 
ponents. Following Refs. p^ , |l8|, a Gilbert damping with damping 
parameter e was chosen in Eq. ( |l.42| ), chiefly because it is isotropic, in 
contrast to the Landau-Lifshitz damping ffStl. 



Strictly speaking, the three equations ( 1.42 ) do not really represent 
Langevin equations, because all the components of dS"'- /dt appear in each 
equation due to the cross product. To properly introduce the noise, one 
first has to group all the time derivatives on the 1. h. s. of the equation, and 
only then one can introduce independent white noise terms for each spin 
component. But this procedure eventually produces only a renormalization 
of the damping parameter e in the order of . In the simulations values of 
e in the order of 10~^ were used, and therefore the correction factor can 
even be neglected. 

Another, even more important issue is the interpretation of the stochastic 
differential equation (1.42). As the noise is multiplicative, Ito and Stratono- 
vic interpretations do not coincide, in contrast to the additive noise case pO| , 
In principle, when thinking of thermal excitations interacting with the 
spins, there would be a finite correlation time which would lead to a colored 
noise term. Taking white noise means taking the limit of zero correlation 
time, and therefore the stochastic Landau-Lifshitz equation (1.42) must be 
interpreted in the Stratonovic sense. 

As noted in the beginning of this subsection, the remaining three steps 
of the stated procedure are straightforward. Therefore we present here only 
the results for the case of non-planar vortices. 



Step 2: Stochastic equation of motion 



The generalized travelling wave ansatz ( 1.29 ) up to order n 
the 3'''^-order stochastic o. d. e. 

AX + MX + GX = F(X) + F™""(t) 



2 yields 



(1.45) 



The tensors A,M and G are the same as in Eqs. ( 1.19D , (1.25) and ( 1.35| ), 
except that all the vanishing components in these three expressions are now 
replaced by damping terms proportional to e. Thus the damping appears 
at every order in Eq. ( 1.45| ). Without F™"'*, the solution is obtained in 
complete analogy to section E^: Two cycloids are superimposed on an 
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outward spiral on a circle (Fig. g). In the simulations the purpose of the 
damping is to dissipate the energy which is supplied to the system by the 
kicks of the noise. Therefore the range of e (for a given system size L) must 
be determined in which the cycloidal frequencies a;i^2 in the trajectories are 
not influenced by the damping p^ . The result is the condition eL <C 6. 
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FIGURE 8. Schematic sketch of the vortex motion as governed by the Lan- 
dau-Lifshitz equation with Gilbert damping. The plot is approximate and does 
not correspond to an actual simulation. 



The stochastic force in (1.45) 

Fr''{t)-^ |d2r^h(r,t), (1.46) 

which stems from the multiplicative noise in Eq. ( |1.42| ), has zero mean and 
the correlation function 

Putting the vortex on the center of a circular system of radius L with free 
BC, the leading contrib ution to the v ariance is obtained |Q by using the 
static vortex structure ( l.lSj ) to (|l.l7 ) 

Yav (Ff) ^ Dv6{t-t') (1.48) 

with the vortex diffusion constant 

Dv = Dtt fin — + C ) , (1.49) 
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here C stems from the inner vortex region < r < and is obtained by 
numerical integration p7|. 



The variance (1.481) is a remarkable result for two reasons: 



(a) The stochastic forces acting on the vortex represent additive white noise 
with an effective diffusion constant. The point is that this was shown by 
going from the microscopic level, where multiplicative white noise was im- 
plemented, to the level of the collective variables. This approach is much 
more satisfying than assuming ad hoc a noise term on the collective variable 
level. 



(b) The mean and variance of the stochastic force ( lAi ) turn out to be the 
same as those of the stochastic force 

^f'' = ^2 /d'^ (s X ^) 77(r, t) , (1.50) 

which resulted from starting with additive white noise J7(r, t) in the Landau- 
Lifshitz equation; this noise has already been discussed in step 1. This 
equivalence of the two forces means that in the case of the vortices additive 
noise and the multiplicative noise in Eq. ((TH) have the same effect on the 



vortex dynamics. This is a nontrivial result, because normally a qualitative 
difference is expected (cf. the discussion in step 1). Qualitatively, this result 
can be understood by noting that the renormalization of the spin length in 
the case of additive noise introduces effectively a multiplicative noise term, 
in addition to the additive one (cf. footnote in step 1). 

Step 3: Variances of the vortex trajectory 

The inhomogeneous stochastic 3'''^-order o. d. e. (1.45) can be solved 



by a standard Green's function formalism, because the 1. h. s. is linear and 
the deterministic force F(X) on the r. h. s. can be linearized by expanding 
in the displacement x(t) = X(i) — X°(t) from the mean trajectory X*'(t) on 



which the vortex is driven by F, cf. Eq. (l.Sf). Knowing x(<), the variance 
matrix 

(T%{t) = {x^Xj) - {x^){xj). (1.51) 

can be calculated; here i, j denote polar coordinates R and 0. Each element 



of the matrix (1.51) turns out to consist of 36 terms. In order to facilitate 
the discussion several simplifications can be made which give the following 
resuhs BtII: 



7 , . D 



(2^)2 



J_ (1 _ e-2/3*) _ 1e-^' sincot + —e'^^' sm2iut 
4/3 Lu 4w 



(1.52) 

with LU ~ y/uJiLU2 and ~ e/5 is the damping constant for the vortex mo- 
tion. For large times, t ^ 1//3, only the first term remains and the variance 
is linear in time, the standard random walk result. Interestingly, this result 
is identical to the one which is obtained by omitting the 2"*^- and 3'''*-order 
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terms in the stochastic equation of motion (1.45). This obviously means 



that these two terms have the foUowing effects: 

(a) They produce the oscillatory parts in Eq. ( 1.52| ), which are naturally 
connected to the cycloidal oscillations in the vortex trajectories, 

(b) for small times, t ^ 1//3, the slope of crlj^, averaged over the oscilla- 
tions, is larger by a factor of 3/2 compared to late times. 

However, the first effect cannot be observed in the simulations because 
the oscillations are hidden in strong discreteness effects (see next step) . For 
this reason we discuss only the long-time behavior of the other elements of 
the variance matrix 



1 Di 



2 (27r)2 2tt 
Dv 



(2^)2 



1 

^+3 



Mo 
27r 



(1.53) 
(1.54) 



with k = l- F„/{F^Ro) and Fq = F{Ro), = F'{R„). The quadratic and 
cubic time dependences are non-standard results which appear in addition 
to the standard linear dependence because the deterministic force field 
F{R) = Fq + Fq{R — Rq) is inhomogeneous. _F is a radial force which drives 



the vortex in the azimuthal direction, due do the gyrocoupling force (1.3) 



Therefore, only the ^-components of a are affected: aj^^ acquires a factor 
fcFQ/(27r) • t, while u^i acquires it twice. 



Step 4- Langevin dynamics simulations 



1.42) was solved numerically 



The stochastic Landau-Lifshitz equation 
for a large circular lattice with free BC and one vortex driven by a radial 
image force ^ , 
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47 1 . The mean trajectory is an outward spiral. Therefore 



a small damping parameter e was chosen in order to allow long integration 
times. 

There are several qualitatively distinct temperature regimes. For < 
T < T3 « 0.05 (in dimensionless units, where Tkt ~ 0.70 for the XY-model 
52 , 53)), two frequencies are observed in the oscillations around the mean 
traject ory which can be identified with the cycloidal frequencies wi^2 in 
section p,.3\ As these frequ encie s are constant in the whole regime, the 3'''^- 
order equation of motion ( 1.45| ) with temperature independent parameters 
can in fact describe the vortex dynamics. 

For Tj, <T <TiK, 0.3 the above two frequencies cannot be observed any 
longer due to large fluctuations , thus the f^'-order equation of motion is 
sufficient here. However, in some of the runs the vortex suddenly changed 
its direction of motion; this will be explained in the next subsection. 

Finally, for T > Ti , a single vortex theory is no longer adequate because 
here the probability for the spontaneous appearance of vortex-antivortex 
pairs in the neighborhood of the single vortex becomes too large. 



1. Dynamics of Vortices in Two-Dimensional Magnets 25 




FIGURE 9. Vortex diffusion constant Dv as a function of ter npera ture, for 
e = 0.002 and L = 24. Solid line: Theoretical results from Eq. ( 1.49 ); dashed 
line: Adjusted Dv from fitting the theoretical curves for a^{t) to the simulation 
data. 



The linear, quadratic, and cubic time dependences in Eqs. ( 1.52 ) - ( 1.54 ) 
are well confirmed by the simulations. Therefore the factor Dy can be 
fitted and turns out to differ from the predicted vortex diffusion constant 
(1.49) only by a nearly temperature independent factor of about 1.8, see 
Fig. y (The constant C in (|l.49D was obtained by a numerical integration 
over the vortex core ]46| ). This agreement is amazingly good, taking into 
account that the simulations were performed for a discrete lattice while the 
theory works in the continuum limit and uses additional approximations, 
like the expansion ( |l.36| ) of the Coulomb force. 




FIGURE 10. Variance of the radial vortex coordinate averaged over 2000 real- 
izations, for T = 0.003, e — 0.002, and L — 24. The dashed vertical lines indicate 
the times at which the vortex center moves over ridges of the Peierls-Nabarro 
periodic potential. 
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Although discreteness effects are hardly visible in the trajectories of non- 
planar vortices at zero temperature see also Fig. ^ these effects are 
very important at finite temperatures: Fig. |o| clearly demonstrates that the 
variance of the trajectory shows a pronounced peak whenever the vortex 
center moves over a ridge of the Peierls-Nabarro potential. Due to this 
effect, the predicted cycloidal oscillations (1-52) in the variance cannot be 
observed. 



3.3 Noise-induced transitions between opposite polarizations 

As already mentioned above, a vortex in a thermal bath can suddenly 
change its direction of motion on its outward spiral. A closer inspection 
shows that this occurs because, opposite to the vorticity q, the polarization 
p is not a constant of motion for a discrete system: The out-of-plane vortex 
structure can flip to the other side of the lattice plane due to the stochastic 
forces acting on the core spins. Then the direction of the gyrovector ( |l.4| ) is 



reversed and according to the Thiele equation (1.6) the direction of motion 
is reversed as well (The same holds for the 3'''*-order equation ( 1.3l| ) because 
A ~ G in Eq. (p|)). 



In a preprint ]54| ] the cores of both planar and non-planar vortices are 
described by a discrete Hamiltonian, similar to the one which was used 
for the study of the instability at S = Sc Using the stochastic 



Landau-Lifshitz equation (1.42) the Fokker-Planck equation is derived. Its 
stationary solution exhibits two maxima for the two possible polarizations 
of the non-planar vortex and a saddle point for the planar vortex, if the 
anisotropy parameter i5 lies in a certain temperature dependent range. The 
rate k for the transition from one polarization to the opposite one is calcu- 



lated in analogy to hanger's instanton theory |5^, 56 1, using the fact that 
for 6 Sc there is a soft mode among the normal modes which were ob- 
tained numerically for a system with one vortex Taking into account 
only the four innermost spins of the core, a very simple result is obtained 

K ^ ^^{1 - S)^ - {1 - Sc)^ (1.55) 

ZTT 

where AE is the energy difference between the planar and the non-planar 
vortex. For S — 0.1, the average transition times r = k^^ are 100440, 
4837, and 1060 for T = 0.1, 0.15, and 0.2 respectively. Despite of the 
crude model for the vortex core, these values agree rather well with the 
transition times from Langevin dynamics simulations: t — 92516, 4016, 
and 825 for the above temperatures, with statistical errors of 22%, 6%, and 
10%, respectively. 
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4 Dynamics above the Kosterlitz-Thouless 
transition 

4.1 Vortex-gas approach 

This is a phenomenological theory which is based on the foUowing assump- 
tions: Above Tkti the free vortices form a dilute gas and move either dif- 
fusively lls^ , [5^ , or ballistically In the former case, only spin auto- 
correlation functions were calculated which lead to dynamic form factors 
without wavevector dependence. Therefore we review only the latter case 

The density n„ and the r. m. s. velocity v are the only free parameters. 
According to Kosterlitz and Thouless ^ 

^ (1.56) 



(20 



where ^ is the static spin correlation length which diverges with an essential 
singularity for T Tkt from above. Therefore the vortex gas is in fact 
dilute if T is close enough to Tkt- 

The vortex density is homogeneous only on the average, locally the dis- 
tribution is expected to be random. Therefore the net force is zero on 
the average and the distribution around zero is Gaussian, which yields a 
Maxwellian velocity distribution. This also holds for the non-planar vor- 
tices because the velocity is proportional to the cross product of the net 



force and the gyrovector, due to the Thiele equation (1.6) 



However, the assumption of a ballistic vortex motion is problematic (sec- 



tion 4.3). Under this assumption the dynamic spin-spin correlation func- 
tions 

So.a{r,t) = {Sc,{r,t)So.{0,0)),a = x,y,z, (1.57) 

and their space-time Fourier transforms, namely the dynamic form factors 
Saai^,^), can be calculated analytically | |5^ . 

There is an important difference between in-plane correlations (a = x,y) 
and out-of-plane ones (a = z): As the in-plane vortex structure is not 
localized (it has no spatial Fourier transform), Sxx = Syy are only globally 
sensitive to the presence of the vortices and the characteristic length scale 
is their average distance 2^. The dynamic form factor exhibits a "cenfraZ 
peaA;" , i. e. a peak around u = 0, with the (squared) Lorentzian form 

Sxx{q,u;)^—- J/^ (1.58) 

where 7 = y/Trv/{2£^). Here the core structure did not enter the result, 
because the theory was worked out on a length scale much larger than 



the core radius ry. The integrated intensity Ix{q) of (1.5?) is inversely 
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proportional to the density , therefore the motion of the vortices actually 
destroys correlations. 

For the out-of-plane correlations Szz the situation is completely different: 
The non-planar vortices have statically the localized Sz structure (1.16), 
( 1.17| ) which has a spatial Fourier transform, namely the static form factor 
/(g). Therefore Szz is locally sensitive to the vortices, i. e., to their size and 
shape. Consequently 



Szzin^i^) 



\m\' 



47r5/2TJ 



q 



exp 



(1.59) 



contains f{q) and is proportional to the density. This Gaussian central peak 
simply reflects the Maxwellian velocity distribution. 

For the planar vortices the situation is more complicated because t hey d o 
not h ave a static Sz structure but only a dynamic one, namely Eqs. ( 1.23 ), 
(1.24). Therefore the form factor is velocity dependent which yields a 
more co mpli cated result for Szz 0, containing the same Gausssian as 
in Eq. (1.59). However, the intensity of this peak is much smaller than 
that of the peak ( 1.59| ), because the dynamic Sz components are about two 
orders of magnitude smaller than the static ones pi] . 



4-2 Comparison with simulations and experiments 

In combined Monte-Carlo (MC) and Spin Dynamics (SD) simulations many 



features of the predicted dynamic form factors (1.58) and (1.5!;) were con- 
firmed: 

(1) In-plane correlations: Both for the XY- model ((5 = 1) ]59| and for the 
weakly anisotropic case {S < Sc) (6^, the observed Sxxi<i, exhibits a cen- 
tral peak for T > Tkt (but not for T < Tkt, as expected). The statistical 
errors were too large to decide about the shape. Howe ver, t he width Tx{q) 
and the integrated intensity Ix{q) of the Lorentzian ( |l.5^ ) can be fitted 
to the data, in this way the free parameters ^ and v are determined. ^(T) 
agrees rather well with the static correlation length . v{T) increases with 
T and then saturates; except of a factor of about 2, this agrees with a result 
of Huber ^ who assumed a diffusive vortex motion and calculated the 
velocity autocorrelation function. 

Central peaks were also observed in inelastic neutron scattering exper- 
iments: For BaCo2(As04)2 [0 and Rb2CrCl4 the measurement was 
performed for only a few g-values. The reported widths of the peaks differ 
from the predicted T^iq) by factors of about 7 and 2.5, resp.. However, in 
this comparison one has to take into account that the theory has neglected 
many features of the real quasi-2D materials: E. g., the lattice structure, a 
pronounced in- plane anisotropy in the case of BaCo2(As04)2, and quan- 
tum effects. 

For C0CI2 graphite intercalation compounds the q-dependence of the 
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central peak width was measured and agrees qualitatively with rj.{q) 
from the peak (1.58). 

(2) Out-of-plane correlations: In MC-SD simulations, Szz{(It^) exhibits a 
central peak only together with a spin wave peak which sits on its shoulder 
(in contrast to S^x where the spin waves are strongly softened for small q) . 

For the XY- model (S — 1), which bears planar vortices, t he in tensity 
of the central peak is expected to be small (see end of section O). Early 
simulations with large statistical errors had difficulties to substract the 
spin wave peak and reported only upper bounds for the width r^(g) and 
the intensity 1^ (q) of the central peak. Recent simulations with much better 
statistics investigated the range Sc < S < 1. Some papers did not find a 
central peak |6^, |5^ others did Q . 

However, for the weakly anisotropic case (6 < Sc) with non-planar vor- 
tices a central peak was indeed found |]63| . Experimentally, S^zio.,^) was 
first not accessible because of intensity problems. Only the use of spin- 
polarized neutron beams made it possible to observe a central peak in 
Rb2CrCl4 and to distinguish it from the spin wave peak and from other 
contributions 1 35|. The measured width is practically equal to the width of 
the Gaussian (1.59). However, such a good agreement seems to be acciden- 
tal because this material exhibits a breaking of the rotational symmetry in 
the XY-plane which is described by a more complicated Hamiltonian | |6^ . 
For the same reason it is not clear whether this material actually belongs 
to the case 6 < Sc- 



4-3 Vortex motion in Monte Carlo simulations 

The vortex gas approach assumes for single free vortices a diffusive |5^, ^ 
or ballistic motion. Very recently this was tested in MC-simulations 
by monitoring the position of each vortex in the system (free or bound 
in a pair) as a function of time ||^, ||, |6|l. The surprising result is that 
a single vortex very seldom moves freely over a larger distance. Normally 
the vortex travels only one or a few lattice spacings until it annihilates 
with the antivortex of a pair which meanwhile appeared spontaneously in 
the neighborhood. Another possibility is that the single vortex docks on a 
pair or a cluster of pairs and after a while one of the vortices leaves the 
cluster. These results confirm suggestions that vortices cannot move freely 
for more than a few lattice spacings, which were made by computing the 
vortex density-density correlation function |]70| , |67| . 

However, the interpretation of these results is not clear at all. One pos- 
sibility is that the vortex gas approach is not valid but then the 
striking qualitative agreement with the central peaks in the simulations 
and experiments is not understood. We favour another possibility: Proba- 
bly only the effective vortex motion is important for the dynamic correla- 
tion function Sxx{'^,t) ( 1.57| ). This would mean that it does not matter if 



a vortex is annihilated with the antivortex of a pair in the considered time 
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interval [0, t], because the vortex of that pair continues the travel instead of 
the original vortex. Only effective lifetimes are seriously affected. This pic- 
ture is strongly supported by the mechanism l59|] which yields the central 



peak Sxxi^lT^) in Eq. (1.58): Looking at the orientation of one particular 



spin in the XY-plane, one sees that its direction is reversed after the vortex 
center has gone over this spin or its neighborhood. Thus Sx and Sy have 
changed their signs after the vortex is gone, and therefore 



(1.60) 



Here N{r,t) is the number of vortices which pass an arbitrary noninter- 
secting contour connecting (0,0) and (r,t). Obviously, it does not matter 
if those vortices are replaced by other vortices during their travel (except 
maybe if the replacement happens just when the contour is passed). The 
evaluation of Eq. (1.60) leads to the central peak (1.58) |59 . 



5 Conclusion 

The vortex dynamics at zero temperature can be well understood by a 
collective variable theory: A generalized travelling wave ansatz which al- 
lows for deformations of the vortex shape due to velocity, acceleration etc., 
leads to equations of motion. In the case of non-planar vortices, the trajec- 
tories exhibit a superposition of cycloids which is fully confirmed by spin 
dynamics simulations. 

This collective variable theory can be generalized to finite temperatures 
which yields stochastic equations of motion. The vortex motion is diffusive 
and agrees well with Langevin Dynamics simulations. Moreover, the rate of 
noise-induced transitions between vortex states with opposite polarization 
is calculated and agrees with the simulations. 

However, above the Kosterlitz-Thouless transition temperature there is 
so far no theory which can fully explain the central peaks which were ob- 
served both in inelastic neutron scattering experiments and in combined 
Monte-Carlo and spin dynamics simulations. A qualitative agreement is 
achieved by a phenomenological vortex gas theory, but one of its assump- 
tions, namely ballistic vortex motion, is questionable. Probably, both the 
diffusive character of the vortex motion and annihilation and creation pro- 
cesses must be incorporated into a theory which can fully explain the above 
facts. It is likely that the situation is similar for all Kosterlitz-Thouless (and 
many other) phase transitions. 
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